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We study the effect of interactions between objects floating at fluid interfaces, for the 
case in which the objects are primarily supported by surface tension. We give conditions 
on the density and size of these objects for equilibrium to be possible and show that two 
objects that float when well-separated may sink as the separation between the objects is 
decreased. Finally, we examine the equilbrium of a raft of strips floating at an interface, 
and find that rafts of sufficiently low density may have infinite spatial extent, but that 
above a critical raft density, all rafts sink if they are sufficiently large. We compare our 
numerical and asymptotic results with some simple table-top experiments, and find good 
quantitative agreement. 



1. Introduction 

A common table-top demonstration of the effects of surface tension is to float a metal 
needle horizontally on water: even though the density of the needle is much greater 
than that of water, the needle is able to remain afloat because of the relatively large 
vertical compon ent of surface ten sion. This effect is a matter of life or death for water- 
walking insects l|Bush fe Hull20of^l . and is also important in practical settings such as 
the self-assembl y of small metallic components in to macroscopic structures via capillary 
notation forces l|Whitesides fc Grzvbowskill2002j) . In this engineering setting an object 
should not only float when isolated at the interface, but must also remain afloat af- 
ter it has come into contact with other interfacial objects, and portions of the menis- 
cus that supported it have been eliminated. Although the interactions that cause in- 
terfacial objects t o come into contact and form clusters have been studied extensively 
(see, for example. iMansfield. Senangi fc Eastwoodlll997t iKralchevskv fc NagvamalfeoOft 
IVella fc MahadevanlbOOIil) . the implications of such interactions on the objects' ability to 
remain afloat have not been considered previously. 

Here we consider the effects of these interactions via a series of model calculations that 
shed light on the physical and mathematical concepts that are at work in such situations. 
For simplicity, the calculations presented here are purely two-dimensional, though the 
same physical ideas apply to three-dimensional problems. 



2. Two horizontal cylinders 

Perhaps the most natural way to characterise the effects of interaction is to ask how 
the maximum vertical load that can be supported by two floating cylinders varies as the 
distance between them is altered. We thus consider two cylinders of infinite length lying 
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Figure 1. Cross-section of two parallel, horizontal cylinders lying at an interface with a 
non-dimensional centre-centre separation of 2A. 

horizontally at the interface between two fluids of densities pa < PBi as shown in figure 
^ We assume that these cylinders are non- wetting so that the contact angle 9, a property 
of the three phases that meet at the contact line, satisfies 9 > tt/2. 

We non-dimensionalise forces per unit length by the surface tension coefficient, ~/ab, 
and lengths by the capillary length, £ c = (^ab/(pb — PAjg) 1 ^ 2 , and use non-dimensional 
variables henceforth. We wish to determine the maximum weight per unit length, W , 
that can be supported by each of two identical cylinders with radius R and centre-centre 
separation 2A. 

To remain afloat each individual cylinder must satisfy a condition of vertical force 
balance: their weight (or other load) must be balanced by the vertical contributions of 
surface tension and the hydrostatic pressure acting on the wetted surface of the cylinder. 
We assume that an external horizontal force is applied to maintain the separation of the 
cylinders and so do not consider the balance of horizontal forces explicitly. 

Using the notation of figure ^ the vertical force balance condition may be written 
W=U 1 + U 2 where 

Ui = - sin(0 + ipi) - H Rwiipi + ^R 2 ^ + sin?/>i cos V>0 (i = 1,2), (2.1) 

are the contributions to the vertical upthrust provided by the deformation on each half 
of the cylinder separately, and Hq is the height of the cylinders' centres above the unde- 
formcd free surface. Physically, the first term on the right hand side of (|2.1J) is the vertical 
component of surface tension, and the second and third terms quantify the resultant of 
hydrostatic pressure acting on the wetted perimeter of the cylinder. The latter is given 
by the weight of water that would fill the dashed area in figurelTI (see lKelleilll998j) . 

The angles ip\ and 4>2 are determined by the interfacial shape, which is governed 
by the balance between hydrostatic pressure and the pressure jump across the inter- 
face associated with interfacial tension. This balance is expressed mathematically by the 
Laplace- Young equation. In two dimensions this is 

H xx =H(l + H 2 x f/ 2 , (2.2) 

where H(X) is the deflection of the interface (again measured positive upwards) from 
the horizontal, and subscripts denote differentiation. Since the exterior meniscus extends 
to infinity, the first integral of l|2.2|) is particularly simple in this instance and allows the 
height of the contact line, , to be related to the interfacial inclination, <j>, via 

cos0= 1 - \Bl. (2.3) 

This, together with the geometrical condition <fi = 9 + tpi — ir, allows ipi to be eliminated 
from (|2.1|) in favour of H (= H* + Rcosipi) and 9. 

For the interior meniscus, we simultaneously obtain tp 2 and the shape H(X), by using 
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Figure 2. Numerical results for the maximum load that can be supported by a single cylinder 
in the presence of another a distance 2 A away when 9 = 2-7r/3 for several values of the Bond 
number, B = R 2 . The dashed line shows the linear approximation 12. 51 for the limiting case 
B = when A < 1. 



the MATLAB routine bvp4c to solve the nonlinear eigenproblem 



on [i?sin ip 2 , A]. 

With the angles ipi and ip2 calculated, W(Hq) can be determined from JZ3J, and the 
maximum load that can be supported, W m ax, can be found numerically by varying Hq. 
Of particular interest is the dependence of lf max on the cylinder separation, which is 
shown for several values of the Bond number B = R 2 in figure [21 This plot includes the 
limiting case B = 0, corresponding to the application of two point forces to the interface. 

The results presented in figure El show that as the distance between two cylinders 
decreases, the maximum vertical load that can be supported by each cylinder decreases. 
Physically, this result is intuitive since even though the interior meniscus is not completely 
eliminated in this instance, the vertical force that this meniscus can exert on the cylinder 
is diminished by the symmetry requirement that Hx{A) = 0. In particular, for small B 
and A the total force that can be supported by each cylinder is around half of that which 
can be supported by an isolated cylinder. This corresponds to the simple physical picture 
that for small B ond number, the restoring force is supplied primarily by the deformation 
of the meniscus l)Hu. Chan fc BusnH 2003): when the interior meniscus is eliminated, the 
contact line length per cylinder, and hence the force that surface tension can provide, 
are halved. From this we expect that very dense objects that float when isolated at an 
interface might sink as they approach one anoth er. Since floating objec ts move towards 
one another due to capillary flotation forces ( see iMansfield et all\l99l[ for example), it 
seems likely that this effect may be ubiquitous for dense objects floating at an interface 
and may also have practical implications. 

For B — we can compute the asymptotic form of W max for A -C 1 by noting that 
for small separations the interior meniscus has small gradients and the Laplace- Young 
equation l|2.2|) may be approximated by Hxx — H, which has the solution H(X) = 
-ffoCOsh(A — A)/ cosh A. Thus, the vertical force provided by the deformation is W = 
- J ff (tanhA+(l-iJ 2 /4) 1/2 ), which is extremised when (H Q 2 ~2)/(4-H 2 ) 1 / 2 = tanhA. 
Choosing the real root of this quartic corresponding to a maximum in W and making 



Hxx 
H x (Rsmip 2 ) 
H(R sin ip 2 ) 
H X (A) 



H{l + H 2 x f /2 : 
tan(6» + ^ 2 ), 
H — Rcos4>2, 
0, 



(2.4) 
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Figure 3. Cross-section of two shallow, touching strips floating at a liquid-fluid interface. 

consistent use of A <C 1, W ma x can be expanded as a series in A. We obtain 

^max = l + V2A + 0(A 2 ), (2.5) 
which compares favourably with the numerically computed results presented in figure El 

3. Two touching strips 

Whilst the scenario considered in the previous section may be relevant in practical 
situations, it does not lend itself to particularly simple experimental validation. To allow 
for such a comparison, we now consider the equilibrium of two infinitely long, shallow 
strips of dimensional thickness £ c t, width 2£ c [3, and density p s , floating with their long 
edges in contact so that the interior meniscus is completely eliminated. The configuration 
is shown schematically in figure |3] Here, we are no longer bound by a contact angle 
condition but instead assume that the meniscus is pinned to the uppermost corners of 
the strips. The additional complication of the strip's angle of inclination to the horizontal, 
a, is determined by the balance of torques. (This condition is satisfied au tomatically for 
shape s with circular cross-section and constant contact angle, as shown bv lSingh &: Heslal 
<l2004lU 

Equating moments about the point of contact (thereby eliminating the need to calcu- 
late the tension force that the strips exert on one another) and balancing vertical forces, 
we obtain the conditions for equilibrium 

T>(3 cosa — sin(</> — a) — (3(H + |/3sina), (3.1) 

T>/3 = \ sin</> — (3cosa(H + /Jsina), (3.2) 

where 

y^^-P^ (3.3) 

PB ~ PA 

is the appropriate ratio of the density of the strips to those of the surrounding fluids. After 
eliminating T> between (|3.1I) and l|3.2|l and using Ij2.3|l with the relation iJ* = iJo+2/3 sin a 
to eliminate 4>, we have a single equation for a given particular values of /3 and Ho- Thus, 
for fixed (3 and a given value of Hq, we may solve for a and deduce the corresponding 
value of T> from l|3.2|l . By varying Hq we are then able to calculate the maximum value 
of T> for which equilibrium is possible, much as before. The numerical results of this 
calculation are presented in figure 0] 

Also shown in figure 0] are experimental results showing points in (/3, T>) parame- 
ter space for which two identical strips remained afloat or sank upon touching. These 
experiments were performed with strips of stainless-steel shim of length 69 mm with 
p s = 7905 kg m -3 and thickness 0.4 or 0.5 mm. These were floated on aqueous solutions 
of 0%, 10% or 25% methanol in air (so that pa <§C Pb), allowing a wide range of values of 
(3 and V to be probed. The strips wer e then allowed to com e into contact naturally via 
the mutually attractive notation force an sfield et ra/]IT99"7V The data are plotted with 
horizontal and vertical error bars. The former indicate the uncertainty in the measure- 
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Figure 4. Numerically computed values of P ma i as the half-width of the strips, j3, is varied 
(solid line) . Experimental results (as described in text) are shown by points x (strips that sink) 
and O (strips that float). 
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Figure 5. Schematic illustration of a flexible raft floating at an interface. 

merit of the strip widths. The latter indicate the uncertainty in the additional vertical 
force contribution of the ends (since the strips are of finite length), which may be shown 
to be equivalent to an uncertainty in the effective value of T>. The agreement between 
our experiments and theory in this instance is very good. 



4. The floating of a flexible raft 

By adding additional strips to a floating pair of strips, a flexible raft is formed. Given 
the analysis of the preceding sections it is natural to expect that as the raft is lengthened 
in this manner, there will come a point where its weight (which scales with its total 
length) exceeds the force that can be supplied by surface tension (which is constant) and 
so the raft should sink. The situation is complicated by the fact that the raft may bow 
in it s midd le, displacing a considerable amount of liquid in this region, as pointed out by 
ISaiJ l|2002t) . We now address the question of whether, for a raft of given weight per unit 
length, there is a maximum raft length before sinking occurs. 

We tackle this problem by treating the raft as a continuum, shown schematically in 
fig. [SJ and formulatin g an equation for the deformation of such a raft. This generalises 
the linear analysis of lMansfield et al\ l)l997|) and allows us to consider situations in which 
interfacial deformations are no longer small, including the existence of a threshold length 
for sinking. 

4.1. Governing equation 

We use a variational approach to determine the shape (X(S),H(S)) of the raft and 
the surrounding meniscus, though the same result may also be obtained by considering 
the force balance on an infinitesimal raft element. The non-dimensional arc-length, 5, is 
measured from the raft's axis of symmetry at S = 0, with the two ends of the raft being 
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at S = ±L. For simplicity, we n e glect the intrinsic bending stiffness of the raft, although 
IVella. Aussillous fc Mahadevanl l|2004h have shown that interfacial rafts do, in general, 
have some resistance to bending. The variational principle states that raft shapes must 
minimise the energy of the system over variations in H(S) and X(S), subject to the 
constraint that X§ + H s — 1. Introducing a Lagrange multiplier X(S) associated with 
this constraint, we find that equilibrium raft shapes extremise 

£ = ^ (x s (H 2 /2-l)+VH X +(l-x) + HS)[(X 2 s + H 2 s ) 1 / 2 -l])dS, (4.1) 

where V was defined in (|3.3|) and 

■«-ft mil 

is the indicator function of the raft. 

The first two terms in the integral (|4.1() correspond to the gravitational energy of the 
displaced fluid and the raft, the third term is the surface energy of the uncovered liquid 
area, and the final term ensures that the constraint Xg + H s = 1 is satisfied. Note that 
a small increase in arc-length such that X§ + Hg > 1 increases the energy of the system 
so that the Lagrange multiplier X(S) may be interpreted physically as the tension in 
the raft/meniscus. That the raft can support a tension at all may seem counterintuitive. 
It is a consequence of the attractive capillary interaction that would exist between two 
infinitesimally separated raft elements. 

Requiring £ to be stationary with respect to variations in H(S) and X(S) yields 
differential equations for X and H . Using the differential form of the constraint, X$Xss + 
HsHss — 0, we may eliminate A to obtain As = \DHs- This may be integrated using the 
boundary term from integration by parts at ±oo, the boundary conditions H(±oo) = 
and Xs(±oo) = 1 as well as the continuity of A at the raft edge, S — ±L, to give 
A = 1 + xD(H — H*), where = H(±L). We now find the raft shape numerically by 
solving the nonlinear eigenproblem 

X s = cos9, H s = S ine, 9 S H + V™9 



1+V{H-H*Y (4.3) 
X(0) = 0, 6>(0) = 0, 6{L) = 2arcsin(iJ»/2), H(L) = H*, 

for X(S), H(S), 9(S) on [0,L], and H*, using the MATLAB routine bvp4c. The results 
of this computation may be verified by calculation of the quantity 

P{9) = \H 2 + [l + V(H - ff*)] cos 9 - 1, (4.4) 

which is conserved and, from the boundary conditions, equal to 0. 

In the limit o f smal l deformations (|4.3|l reduces to the simpler linear form studied by 
iMansfield et al\ l|l997f) in the context of determining typical raft profiles. Here, however, 
we wish to determine whether a maximum raft length, 2L max , exists and if so find its value 
for a raft of given density T>. To investigate this, small deformation theory is inadequate 
since sinking is an essentially non-linear phenomenon. 

The symmetry condition 9(0) = ensures that > —U/2 and that Hq = H(0) > —T>, 
so that the centre of the raft may sink at most to its neutral buoyancy level. In what 
follows, it will be convenient to treat Hq and V as parameters giving rise to a particular 
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raft semi-length L(Ho,T>); we find 
HZ f 1 

22? ./n fr„ „n, .,,2 r_ „n2l!/2 



{ [2 + HHv - 1)] 2 - [2 - (ff„ + yHZftVy] 2 } 



which follows by changing integration variables from S to i? in L = J L dS*. This allows 
us to consider the behaviour of L for a given value of T> as Hq is varied. 

The tension at the midpoint of the raft is given by 1 — Hq/2, showing that the raft 
goes into compression if Ho < —a/2. Physically this is unrealistic, corresponding to a 
divergence in 6s- li T> < a/2, this situation is avoided automatically since Ho > —T> > 
—a/2 but for V > v2 we must consider this possibility; we therefore consider these two 
cases separately. 

4.2. The case V < a/2 

When V < a/2, the centre of the raft may reach its neutral buoyancy depth Hq = —T> 
without going into compression. Numerical computation of the integral l|4.5|l suggests that 
rafts grow arbitrarily long as Hq \ —T> (see figureEt). To show that this is the case, we 
consider the asymptotic behaviour of the integral i|4.5[l in the limit e = T> + Hq <C If. 

s gions [0, 5] a r 
; lHiricblll99fl . 

these two regions, the two integrands may be simplified using approximations compatible 
with this gearing of S, and the resulting integrals evaluated analytically. Upon expanding 
these results for 5 -C 1 the leading order terms in S cancel, yielding 

L = -/A log (^r) + Mlog f 



This is done by splitting the range of integration into two sub-re gions [0, S] a nd (5,1], 
where 6 is unspecified save for the condition that e < <5 < 1 fsee lHinchlll99C| . Within 



V) ' b V\/2a*(7 + ^ 2 ) 1/2 + 4 




where fj, = (1 - X> 2 / 2 ) 1/2 - This result compares favourably with the numerical results 
in figure HJa). In particular, notice that L diverges logarithmically as Ho \ ~T> (i.e. as 
e — > 0) so that rafts of arbitrary length are possible. It also interesting to note that l|4.6|l 
may be inverted to give an estimate of Ho = —T> + e for given values of D and L — a 
useful result when calculating raft shapes for large L. 

That a raft of sufficiently low density can grow arbitrarily large in horizontal extent 
without sinking seems surprising at first glance. However, as new material is added to 
the raft, it may be accommodated at its neutral buoyancy level without the raft going 
into compression. Therefore, the raft's ability to remain afloat is not jeopardised and it 
is almost obvious that these low density rafts may grow arbitrarily long without sinking. 

4.3. The case V > a/2 

In this case, the raft cannot reach its neutral buoyancy level, invalidating the argument 
just given to explain why, with T> < a/2, rafts may be arbitrarily large. We thus expect 
that a maximum raft length does exist and, further, that the limiting raft has Ho = —a/2. 
Numerical computation of L as a function of Ho indicates that a critical half-length L max 
does exist, but that it is not attained with exactly this value of Ho- Instead, there is a 
competition between the raft sinking deep into the liquid (to support its weight by 
increased hydrostatic pressure) and having its ends a large distance apart (i.e. lower 

f Note that e > 0, since H > -V. 
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Figure 6. (a) Numerical results of the calculation of L as a function of e = Ho + T> (solid 
line) compared to the asymptotic result 14.61 for e < 1 (dashed line) for the case T> = 1.4. 
(6) Main figure: Numerical results of the calculation of L max as a function of the density ratio 
T> > \/2 (solid line), together with the large D asymptotic result Z/ max ~ 1/2? (dashed line). 
Inset: Rescaled graph comparing the numerically computed values (points) of L max with the 
asymptotic expansion 14.91 (solid line). 



pressure but over larger horizontal distances), and some compromise is reached. Given 
the abrupt change in behaviour observed as T> increases past y/2, we are particularly 
interested in the nature of this transition. Numerical computations suggest that for rj 2 = 
V- V2< 1, L max occurs when Hq = —\/2 + crj 2 for some constant c. Motivated by this 
observation, we let Ho = — V2 + crj 2 and again split the domain of integration in l|4.5|l 
into two regions [0,6'} and [S',1] where i) 2 C f « 1. This allows us to calculate L to 
leading order in 77, yielding 



2v / 2arctan —— 



V7\ 2 3 / 4 c 



3 / ' ''(c+l) 1 ^ 



0(77 2 ), 



,2(c+l)y c V 2 ( c + 1 ) 

(4.7) 

where K(jfe) = J* /2 (l - k 2 sin 2 ^)~ 1 / 2 d(j) and E(fc) = f* /2 (l - k 2 sin 2 0) 1 / 2 d0 are the 
complete elliptic integrals of the first and second kinds, respectively. The coefficient of 7/ 
in (|4.7(l has a maximum for fixed 77 at c = c*, where c* satisfies 



K l2^J- 2 H^J- (48) 

Hence c* ~ 0.5332, and we obtain the asymptotic expression 

1 /2 

L max = 2\/2 arctan f ^ - 3.1525 fx> - v^J " + O (v - Vfj, (4.9) 

which compares very favourably with the numerically computed values of L ma x presented 
in the inset of figure HJfr)- 

For the limiting case T> = v2, the above analysis breaks down since then r\ = and 
we lose the freedom to vary Ho. However, by letting e = crj 2 (so that Ho = e — v2) we 
take the limit rj — > of (|4.7(l with e <C 1 fixed to find 

L(e) = 2V2arctan (^) + e 1 /2 2 3/4 |" K ^j _ 2E(^)] + O(e). (4.10) 

This has a maximum value of 2-\/2arctan(v7/3) at e = 0, which is the same value as 
that found from (|4.9(l in the limit V \ It is also reassuring to note that, as T> /* \[2 
with e fixed, the expression in H4.6(l also gives L = arctan(\/7/3) + 0(e 1 ^ 2 ). 

For completeness, we consider finally the limit T> ^> 1. To leading order in I? -1 , the 
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Figure 7. Comparison between experimental realisation of a two-dimensional raft (viewed 
through the side of a transparent tank) and the theoretically predicted shape (superimposed 
white line). The rafts float at an air-water interface and have varying values of T> and L: (a) a 
complete raft with V = 1.02, L = 4.03 (6) one half of a raft with V = 1.27, L = 1.47 and (c) one 
half of a raft with T> = 1.27, L = 2.57. The typical width of each individual strip is 2 mm. The 
black region apparently above the raft is in fact a reflection of the black base of the confining 
tank from the meniscus at the edge of the tank 



integral for L(Ho,T>) is given by 

This has a maximum value of at Ha — —y/2 so that in the limit 2? 3> 1, L max ~ 2? _1 . 
This is precisely as we should expect physically since large density objects can only float 
when the contribution of surface tension dominates that of the buoyancy due to excluded 
volume and, in particular, it must balance the weight of the raft. This asymptotic result 
compares favourably with the the numerical results presented in figure . 

4.4. Comparison with experiment 

A direct comparison between the theoretical results outlined so far and experimental 
results is difficult since we have modelled the raft as a perfectly flexible continuum body 
of infinite extent along its axis of symmetry. Despite these limitations, the theoretical 
raft shapes calculated via this model are in good agreement with simple experiments in 
which thin strips of stainless steel shim are laid side-by-side at an air-water interface, as 
shown in figure — even when the raft consists of only a small number of strips and we 
might not expect the continuum approximation to be valid. 

Although this agreement is encouraging, our main interest lies more in whether there is 
a maximum length for such a raft to remain afloat, as predicted by the model. Practical 
considerations mean it is difficult to produce strips of stainless steel shim narrower than 
about 2 mm in the workshop, so the comparisons we are able to draw between our model 
and experiments can only be semi-quantitative. In spite of these limitations, we find that 
for stainless steel strips of length 69 mm and thickness 0.5 mm the maximum raft-length 
is 4 — 6 mm for an aqueous solution of 25% methanol (so that 1.645 >T>> 1.580) and 
6 — 8 mm for 15% methanol (so that 1.494 > V > 1.424), which are certainly consistent 
with the corresponding theoretical results of 4.6 mm < L max < 4.8 mm and 6.5 mm 
< imra < 7.2 mm, respectively. Here the length was increased by floating additional 
strips near the raft and allowing them to come into contact via the mutually attractive 
capillary flotation forces until the raft was no longer stable and sank. With T> = 1.02 and 
T> = 1.27, we were able to add many strips without any sign of the raft sinking indicating 
that this process might be continued indefinitely. 
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5. Discussion 

In this article, we have quantified the conditions under which objects can remain 
trapped at a fluid-fluid interface, and shown that when the deformation of the meniscus 
is suppressed by the presence of other objects the supporting force that can be generated 
decreases dramatically. For two small, parallel cylinders or strips, the maximum force 
that can be supported close to contact is only that provided by the contribution from the 
exterior meniscus and so sufficiently dense objects sink upon contact. A two-dimensional 
raft of touching, floating strips may compensate partially for this loss of meniscus by 
sinking lower into the fluid. For V < a/2, this effect allows rafts of arbitrary length to 
remain afloat. For T> > \/2, there is a maximum length (dependent on T>) above which 
equilibrium is not possible. 

Although the agreement between the experiments and theory presented here is good, 
our analysis was confined to two dimensions, whereas experiments must be carried out 
in the three-dimensional world. Similarly, we have limited ourselves to considering the 
equilibrium of objects at an interface. We are currently studying the dynamics of sink- 
ing for the case of two touching strips considered in section |3 and find that a simple 
hydrodynamic model produces good agreement with experiments. 

We are grateful to David Page-Croft for his help in the laboratory and Herbert Huppert 
for comments on an earlier draft. DV and RJW are supported by the EPSRC. PDM 
gratefully acknowledges the financial support of Emmanuel College, Cambridge. 
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